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Abstract 

We generalize Huberman-Rudnick universal scaling law for all periodic windows of the logis- 
tic map and show the robustness of (/-Gaussian probability distributions in the vicinity of chaos 
threshold. Our scaling relation is universal for the self-similar windows of the map which exhibit 
period-doubling subharmonic bifurcations. Using this generalized scaling argument, for all periodic 
windows, as chaos threshold is approached, a developing convergence to g-Gaussian is numerically 
obtained both in the central regions and tails of the probability distributions of sums of iterates. 
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I. INTRODUCTION 



It is well-known that many complex systems exhibit transitions from periodic motion 
to chaos through period doubling route like Rayleigh-Benard system in a box forced 
pendulum logistic map [3], Chirikov map [4J etc. The logistic map, defined as 

x t+1 = 1 - ax 2 t , (1) 

(where < a < 2 is the control parameter and the phase space x t is between [—1,1] 
with t — 0, 1, 2, ...), is a good example to observe Feigenbaum route generated by pitchfork 
bifurcations and its universal features. This map has its critical point, denoted by a c , at 
a c = 1.401155189..., which can be approached from left (i.e., from periodic region) via period 
doubling procedure where 2°° periods accumulate at this critical point usually described as 
chaos threshold. This point can also be approached from right (i.e., from chaotic region) 
via band merging procedure where infinite number of bands merge at the critical point jsj]. 
A sketchy view of these approaches from left and from right to a c is given in Fig. UK In 
the chaotic region (a > a c ), there exist many windows of higher periodic cascades of s2 n 
periods, where s is an integer and n is the degree of period-doublings. It is worth noting 
that, within each window, reverse bifurcations of s2 n+l bands merging into s2 n bands can 
easily be detected although the width of these windows decreases rapidly as long as s ^ 1. 
As an example, s = 3 case (namely, period 3 case) is illustrated in Fig. [Tb. 

This self-similar structure of the map is of sensible importance and leads Feigenbaum 
to develop a scaling theory for the non-chaotic period-doubling region of the bifurcations, 
which enables one to localize the control parameter value at the bifurcation from 2 n period 
to 2 n+1 via scaling relation 

|a-a c |~<T r \ (2) 

where 5 = 4.699... is the Feigenbaum constant (6j. On the other hand, it is possible to obtain 
Eq. ([2]) from the famous Huberman-Rudnick scaling law [7|], which shows that the envelope 
of the Lyapunov exponents near a c exhibits a universal behavior, similar to that of an order 
parameter close to the critical point of phase transition. This relation can be written as 

A = Ao|a — a c \ u , (3) 
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FIG. 1: (a) Sketchy view of the standard bifurcation diagram of the logistic map. Approaching 
chaos threshold (a c = 1.401155189...) from left via period-doubling route is evident with bifurcation 
points denoted as a n (n = 1, 2, oo). Corresponding band merging approach to a c is depicted with 
merging points denoted as a n at each of which 2™ bands appear, (b) Sketchy view of period 3 win- 
dow in the chaotic region of the logistic map. Approaching chaos threshold (a c = 1.779818075...) 
of this periodic window from left via period-doubling route is evident with bifurcation points de- 

(3) 

noted as a n . Corresponding band merging approach to a c is depicted with merging points denoted 

aS dyi 

at each of which 3x2" bands appear. 
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where a > a c , v = In 2/ In S, A is the Lyapunov exponent and Ao is a constant. Eq. (J2J) can 
easily be obtained from Eq. ([3]). For a values slightly above the chaos threshold, there exist 
2™ (n = 1,2, ....,oo) chaotic bands, which approach the Feigenbaum attractor as n — > oo 
by the band splitting procedure. In that region, if we start from two trajectories separated 
by a distance <i within one of these chaotic bands, the separation of trajectories increases 
exponentially. If the trajectories start off in one band, after 2 n iterations they will be back 
in the original band. Then they will be separated by the amount 

d 2n = d e X2n = d e x °, (4) 

where Ao = A2™ is the effective Lyapunov exponent (a constant value) Q]. Substituting the 
effective Lyapunov exponent into Eq. ([3]) immediately gives 



2 -n 



\a — a, 



I In 2/ In <5 



(5) 



from where Eq. ([2]) is easily obtained. 

This scaling relation, in fact, is exactly the one used in |8[, where the probability distri- 
butions of the sums of the iterates of the logistic map, as a c is approached from the band 
merging region, have been shown to be well approached by g-Gaussians provided that the 
appropriate number of iterations (N*) is obtained from the above-mentioned scaling relation. 

g-Gaussians, defined as, 



P(y) 



(6) 



P(0) [1 - /3(1 - q)y 2 ]~* for 0(1 - q)y 2 < 1 
otherwise 
(where q < 3 and /3 > are parameters and the latter controls the width of the distribution) 
are the distributions that optimize, under appropriate constraints, the nonadditive entropy 
S g (defined to be S q = (1 — ^i^D / (Q ~~ -0)' on which nonextensive statistical mechanics is 
based (9), [loj|. As q — > 1, g-Gaussians recover the Gaussian distribution. 

Although in Nature many stochastic processes, consist of sum of many independent or 
nearly independent variables, are known to converge Gaussian distribution due to the stan- 
dard central limit theorem fill . [lJ ]. in recent years several complex systems such as low- 
dimensional dissipative maps in the vicinity of chaos threshold 13|4l5|, high dimensional 



dissipative systems 



16] and conservative maps 



3, Q 



are shown to exhibit probability 



distributions that are well approached by g-Gaussians. 



Our main aim in this paper is two-folded: firstly, we try to generalize the Huberman- 
Rudnick scaling law to all periodic windows of the logistic map, secondly, using this gen- 
eralized version of the Huberman-Rudnick scaling law, we analyse the robustness of the 
probability distribution of the sums of iterates of the logistic map as chaos threshold is 
approached. 

II. GENERALIZATION OF THE HUBERMAN-RUDNICK SCALING LAW 

In order to find generalized version of the Huberman-Rudnik scaling law, let us start 

(s) 

by denoting the accumulation point of a particular periodic cycle s as a c . For example, 
and stand for the accumulation points of period 3 and period 5 windows inside the 
chaotic region. These points can be found clS Ojc ^ — 1.779818075... and Oc ] = 1.631019835... . 
Then, one need to check whether the form of the Huberman-Rudnick scaling law is valid for 
all other periodic cycles in the chaotic region. More precisely, one need to check whether 
the exponent v in the scaling law is equal to ln2/ln<5 as in the standard case. Therefore, 
we have first checked this and found that, for all periodic windows, the envelope of the 
Lyapunov exponents, in log A vs log(a — a c ) plot, is given by a slope 0.449, which is nothing 
but In 2/ In 5. Hence we can now write the Huberman-Rudnick scaling law for all periodic 
cycles as 

\ = Xo[a-a^] ln2/hlS . (7) 

At this point, it should be recalled that, for a particular period s window, a trajectory that 
starts in one band will be back in the same band after s2 n (n = 1,2, ...,oo) iterations. If 
we use this feature in the definition of the effective Lyapunov exponent, namely A = As2 n , 
then one can write the Huberman-Rudnick scaling law for other periodic windows as 

\ = \s2 n [a-a^] ln2/lnS (8) 

and 

2- n = s\a-a^\ hl2/hlS . (9) 
This equation enables us to obtain the generalized Huberman-Rudnick scaling law as 
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|a-aW| =r n -&. (10) 

This new scaling relation is valid for all periodic windows including the standard case for 
s — 1, which immediately recovers the standard scaling law given in Eq. (J2J). 



III. ROBUSTNESS OF PROBABILITY DISTRIBUTIONS 

Now let us concentrate on the probability distributions of the sums of iterates of the 
logistic map, which can be written as 

N 

y := {xi ~ (x)) , (11) 
i=i 

where Xj are the iterates of the logistic map and X\ is the initial value regarded as a random 
variable. It has analytically been proved that, for strongly chaotic systems, the probability 



distribution of y becomes Gaussian for iV — > oo [19|, |20j]. Here, the average (...) is calculated 
as time average. On the other hand, as mentioned before in Section I, several complex 
systems of the type low and high dimensional dissipative and conservative exist where the 
probability distribution does not approach to Gaussian, and therefore violating the standard 
central limit theorem due to possible lack of ergodicity and mixing properties. For such 
systems, it is necessary to take the average over not only a large number of N iterations but 
also a large number of M randomly chosen initial values, namely, 

M = ^£X>. e> . (i2) 

3=1 i=l 

As chaos threshold is approached, the logistic map has already been studied in this respect 



21 



221 ] . As it has already been argued in |8|, in principle, in order to attain chaos 
threshold point exactly (i.e., approaching this point with infinite precision), one needs to take 
n — > oo, which, in other words, means that the necessary number of iterations to achieve the 
limit distribution at the chaos threshold is N* — > oo since N* = 2 2n . Since this is, no doubt, 
unattainable in any numerical experiment, one can only approach to this critical point using 
the appropriate values for (a,N*) pairs coming from the Huberman-Rudnick scaling law. 
As long as this scaling law is obeyed, developing g-Gaussian shape of the limit distribution, 
as chaos threshold is approached, has been clearly shown in jg]. On the way of approaching 



chaos threshold, for any approximation level of finite (a, N*) pairs, if the number of iterations 
used is too much larger than iV* and therefore violating the Huberman-Rudnick scaling law 
(i.e., N >> N*), it is of course not surprising that the probability distribution starts to 
approach to Gaussian form from its central part since the system starts to feel that it is 
not exactly at the chaos threshold. Such numerical examples can be found in {g[ 2l|. On 
the other hand, if the number of iterations used is too much smaller than N* and therefore 
violating again the Huberman-Rudnick scaling law (i.e., iV << N*), then the summation 
starts to be inadequate to approach the shape of the limit probability distribution and it 
exhibits a kind of pea ked or multifractal distribution. Such numerical examples have already 



been given in 
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|22|. 

In the remainder of this work, we try to provide further evidence on the robustness of the 
g-Gaussian probability distributions seen as the chaos threshold is approached. In order to 
accomplish this task, we investigate other periodic windows (numerically chosen examples 
are period 3 and 5) making use of our generalized Huberman-Rudnick law. 

As separated band structure for periodic cycle 2 goes from 2 1 to 2°° with 2 n (n = 
1,2, ...,oo), for any periodic cycle s, the same behavior would be to go from s2° to s2°° 
with s2 k {k = 0, 1, oo). It is evident that there is k — > (n — 1) transformations between k 
and n mathematically. Generically, for any periodic cycle s in the chaotic region, one needs 
to perform s 2 k iterations of the map for a given initial value with a control parameter a 
obtained from the generalized scaling law. After s 2 k iterations, the system will basically fall 
into the same band of the band splitting structure. This means that the sum of the iterates 
Ylt=i x i wm essentially approach to a fixed value w = s2 k (x) plus a small correction Aw\ 
which describes the small fluctuations of the position of the s 2 k th iterate within the chaotic 
band. Hence, one can write 

s2 k 

i=i 

If we continue to iterate for another s 2 k times, we obtain 

2s2 k 

y 2 = - (x)) = Aw 2 . (14) 

i=s2 k + l 

The new fluctuation Au> 2 is not expected to be independent from the old one Awi, since 
correlations of iterates decay very slowly if we are close to the critical point. Continuing 



this 2 k times, we finally obtain 



s2 2k 

y 2k = ^2 (Xi- (x)) = Aw 2 k (15) 

i=s2 2k -s2 k +l 

if we iterate the map s 2 2k times in total. The total sum of iterates 

s2 2k 2 k 

y = X)(x i -(x» = X;A«; i (16) 

i=i j=i 

can thus be regarded as a sum of 2 k random variables Awj, each being influenced by the 
structure of the s2 k chaotic bands at distance a — = <5~™~^ from the Feigenbaum 
attractor. At this distance to chaos threshold, in order to see the limit distribution, appro- 
priate number of iterations would be iV* = s2 2k , which corresponds to iV* = s2 2n ~ 2 after 
k — > (n — 1) transformation. 

Now we are ready to check the shape of the probability distribution of any periodic win- 
dows obeying our generalized Huberman-Rudnick scaling law. Chosen examples of possible 
periodic windows are period 3 and period 5 since they are the largest two periodic win- 
dows available in the chaotic region. Although conceptually nothing is changed for small 
sized windows, numerical analysis is getting more difficult as windows sizes are decreasing. 
Numerically used values are given in the Table for our period 3 and 5 analysis. Control 
parameter values a are chosen so that the precision of corresponding n values, coming from 
the generalized Huberman-Rudnick scaling law as 



In 

n = 



a — a c 



(17) 



ln<5 ln2 ' 

would be the same (see Table). This means that we are approaching the critical point with a 
values located on a straight line with a given slope. Our results are given in Fig. 2 for period 3 
and in Fig. 3 for period 5. In both cases four representative points systematically approaching 
the chaos threshold is given. It is clear from Fig. 2a and Fig. 3a that the probability 
distributions of both periodic windows approach to a g-Gaussian. It is also evident that, as 
the chaos threshold is better approached, the tails of the distribution develops better on the 
g-Gaussian, signaling that the limit distribution obtained at the exact chaos threshold point 
would be a g-Gaussian with infinitely long tails. We also present the same data in Fig. 2b 
and Fig. 3b in a different way so that a straight line would be expected for g-Gaussians. 
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l-(l-q)p[yP(0)] 2 

FIG. 2: (a) Probability distribution functions of period 3 (s = 3) for four representative cases 
with different n values. As n increases, numerical convergence to a g-Gaussian with q = 1.63 and 
j3 = 6.3 is appreciated, (b) The same distribution plotted against 1 — (1 — q)(3[yP(0)] 2 on a log-log 
plot for the case which is the closest to the chaos threshold. A straight line is expected with a 
slope 1/(1 — q) if the curve is a (/-Gaussian. It is clearly seen how the straight line is surrounded 
by the log-periodically modulated curves. 
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FIG. 3: (a) Probability distribution functions of period 5 (s = 5) for four representative cases 
with different n values. As n increases, numerical convergence to a g-Gaussian with q = 1.62 and 
j3 = 6.3 is appreciated, (b) The same distribution plotted against 1 — (1 — q)(3[yP(0)] 2 on a log-log 
plot for the case which is the closest to the chaos threshold. A straight line is expected with a 
slope 1/(1 — q) if the curve is a (/-Gaussian. It is clearly seen how the straight line is surrounded 
by the log-periodically modulated curves. 
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Only the closest cases to the chaos threshold for each periodic window are plotted. It is seen 
from these plots that the curves develop on top of a straight line surrounded by log-periodic 
modulations. 



TABLE I: Parameter values used in this work. The values of n obtained from the generalized 
scaling law using Eq. (|17p . the corresponding N* values and the values of q and (5 (estimated from 



simulations) are listed, s = 1 case, already discussed in 
comparison. 



has also been included in the Table for 



s 


a 


2n 


N* 


1 


P 




L40f588 


10.05 


2 io 






ft 


L40f248 


12.05 


2 12 


1.70 


6.2 




1.401175 


14.05 


2 14 








1.40115945 


16.05 


2 16 








1.779819805038384 


14.05 


3 x 2 12 






3 


1.779818446177396 


16.05 


3 x 2 14 


1.64 


6.2 




1.779818155150985 


18.05 


3 x 2 16 








1.779818092822039 


20.05 


3x2 18 








1.63102039110464 


14.05 


5 x 2 12 






5 


1.63101995463619 


16.05 


5 x 2 14 


1.62 


6.2 




1.63101986115802 


18.05 


5 x 2 16 








1.63101984113785 


20.05 


5x2 18 







t Values related to s = 1 case are taken directly from the Table given in ref. [8|. 



IV. CONCLUSIONS 

Our main results obtained in this paper can be summarized as follows: (i) For the logistic 
map having self-similar structure, Huberman-Rudnick universal scaling law has been gen- 
eralized, which becomes now consistent to all periodic windows in the chaotic region of the 
map. This new generalized scaling law is of sensible importance since it enables us to produce 
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self- similar structure of the map and to explain all band merging structures in all available 
periodic windows using only one generalized formula, (ii) The standard Huberman-Rudnick 
scaling law has already been used in js, [23] and g-Gaussian probability distributions have 
been observed as the standard period 2 accumulation point is approached. However, in order 
to test the robustness of g-Gaussian distributions, a first straightforward attempt should be 
to analyse other critical points (chaos thresholds) of different periodic windows located in 
the chaotic region of the logistic map. Since the generalized Huberman-Rudnick scaling law 
obtained in the first part of this paper now enables us to localize appropriate (n,N*) pairs 
as the accumulation point is approached, we managed to check two representative periodic 
windows. For each case studied here (and possibly for all other periodic windows) it is nu- 
merically shown that the g-Gaussian probability distributions with log-periodic oscillations 
are again the observed distributions and developing better as the critical point becomes 
closer. These results clearly indicate the robustness of the g-Gaussian probability distribu- 
tions seen in the vicinity of chaos threshold. Although the obtained q values seem to exhibit 
a slow decreasing tendency as the size of the periodic window becomes smaller, we believe 
that the genuine limit distribution of the chaos threshold (for all s values) would converge 
to a g-Gaussian with a unique q value, which is expected to be in the interval [1.6, 1.75]. 

Finally it is worth mentioning that the results obtained here are expected to be valid 
for all other dissipative maps sharing the same universality class with the logistic map. As 
an open question that can be addressed in a future work, one can mention the analysis of 
appropriate scaling law for the systems exhibiting quasi-periodic route to chaos. 
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